arXiv:l506.02027vl [quant-ph] 5 Jun 2015 


Gauge invariance in simple mechanical systems 

J. Fernando Barbero G., 1,2 ’^ Jorge Prieto, 3, [j] and Eduardo J. S. Villasenor 3,2 ’^] 

1 Instituto de Estructura de la Materia, 

CSIC, Serrano 123, 28006 Madrid, Spain 
2 Grupo de Teonas de Campos y Fisica Estadistica, 

Instituto Universitario Gregorio Millan Barbany, 

Universidad Carlos III de Madrid, Unidad Asociada al IEM-CSIC. 

3 Instituto Gregorio Millan, Grupo de Modelizacion y Simulacion Numerica, 
Universidad Carlos III de Madrid , Avda. de la Universidad 30, 28911 Leganes, Spain 

(Dated: January 9, 2015) 

This article discusses and explains the Hamiltonian formulation for a class of 
simple gauge invariant mechanical systems consisting of point masses and idealized 
rods. The study of these models may be helpful to advanced undergraduate or 
graduate students in theoretical physics to understand, in a familiar context, some 
concepts relevant to the study of classical and quantum field theories. We use a 
geometric approach to derive the Hamiltonian formulation for the model considered 
in the paper: four equal masses connected by six ideal rods. We obtain and discuss 
the meaning of several important elements, in particular, the constraints and the 
Hamiltonian vector fields that define the dynamics of the system, the constraint 
manifold, gauge symmetries, gauge orbits, gauge fixing, and the reduced phase space. 


* fbarbero@iem.cfmac.csic.es 
I jorgeprietoarranz@gmail.com 
1 ejsanche@math.uc3m.es 



2 


I. INTRODUCTION 


Gauge field theories play a central role in the description of the fundamental interactions 
of physics. A popular way to present the concept of gauge invariance is based on the idea of 
turning global symmetries into local ones, involving arbitrary functions, through the intro¬ 
duction of the so called gauge fields. In many contexts gauge theories are defined, more or 
less explicitly, precisely as those obtained by following this procedure. Their physical use¬ 
fulness hinges upon the possibility of hireling observables that are insensitive to the presence 
of these arbitrary elements; i.e. the identification of suitable gauge invariant functions of 
the dynamical variables. 

An indirect consequence of the introduction of local symmetries is the fact that the held 
equations become singular. This singularity manifests itself as the impossibility to solve for 
some of the second order time derivatives of the helds in terms of the other objects present 
in the equations of motion. An associated effect is the possible appearance of arbitrary 
functions in their solutions (notice, though, the existence of singular systems such as the 
Proca held for which no arbitrariness shows up). From this perspective gauge theories are a 
particular instance of the more general models described by singular Lagrangians (i.e. those 
leading to singular Euler-Lagrange equations). 

The traditional way to deal with singular Lagrangians and the canonical quantization 
of the physical models defined by them relies on the ideas developed by Dirac [Tj. A key 
feature of the algorithm introduced by him to get the Hamiltonian formulation for these 
systems is the appearance of constraints, i.e. conditions that the canonical variables must 
satisfy at all times during the evolution. The quantized version of the phase space functions 
that represent these constraints is a key element in Dirac’s approach. 

A common and widespread misconception is to think that gauge theories must necessar¬ 
ily involve helds, reparametrizations and changes of coordinates or reference frames. The 
purpose of this paper is to show that the dynamics of simple mechanical systems, consisting 
of a finite number of point particles connected by ideal rods, can display gauge behavior. In 
order to make our presentation as pedagogical as possible we will focus on a specific example, 
consisting of four equal masses connected by six rods, and compare it with the quintessen¬ 
tial gauge theory: electromagnetism (EM). We will analyze in detail the Lagrangian and 
Hamiltonian formulations for the particular model considered here. As we will show our ex¬ 
ample mimics some of the crucial features of electromagnetism and is richer in some sense, 
in particular regarding its Hamiltonian formulation. 

The best way to understand the essence of Dirac’s construction is in geometric terms. A 
very clear perspective on this issue was provided by Gotay, Nester and Hinds (GNH) [2J 13] 
so we will use their method. Instead of giving an abstract description of the GNH algorithm 
we will introduce it as we perform the actual computations for our model. In our opinion 
the present paper will serve the dual purpose of clarifying some of the concepts behind gauge 
systems in a very simple setting (constraints, gauge orbits, gauge symmetry, gauge fixing, 
reduced phase space...) and also provide a pedagogical introduction to the Hamiltonian 
description of singular systems. 

The structure of the paper is the following. After this introduction we will start to study 
in section [II] a particular, but representative, model that displays the gauge behavior that we 
want to discuss: four particles connected by six rods. Section [Hi] will be devoted to obtaining 
the Hamiltonian formulation for this system in a neat way by using geometric methods 


inspired in the GNH algorithm. The paper ends in section IV with a short discussion 
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A comment is in order here; in order to make the paper accessible to advanced under¬ 
graduate students we are not assuming any prior knowledge of differential geometry on the 
part of the reader (only standard multivariate calculus), however, we will mention by name 
some of the relevant geometric objects to justify their use and show their logical connection 
with the concepts discussed in the paper. We will gloss over several technical points that 
can be skipped in a first approach to this subject. Readers interested in the geometrization 
of classical dynamics are referred to the comprehensive treaty by Abraham and Marsden |3j. 


II. THE FOUR PARTICLE MODEL 

Singular finite dimensional dynamical systems have been considered in some detail in the 
literature (models displaying different pathological behaviors can be found, for instance, in 
the book by Henneaux and Teitelboim [5]). The main drawback of the usual examples is 
their rather artificial character, i.e. they do not describe systems with a simple physical 
interpretation. One of the goals of this paper is to provide one such example. We base our 
approach on the possibility of considering some constraint forces as dynamical variables on 
a par with the standard (generalized) coordinates. The other central idea is to reproduce, 
to some extent, the behaviour of hyperstatic system^] in a dynamical situation. 

The simplest model that we could discuss would be a system of two point masses con¬ 
nected by two ideal rods (massless and completely rigid). It is obvious that the force exerted 
by each rod is undetermined: only their sum can have a physical meaning. Similar arrange¬ 
ments with an arbitrary number N of rods or other collinear models display the same kind 
of behavior. We will discard them here for two main reasons: they are somehow trivial 
on one hand (i.e. the Wrod model with two masses) or non-generic in a concrete sense 
(collinear systems are infinitesimally flexible in the parlance of reference [6]). Furthermore, 
the natural representation of our model, inspired in graph theory, can be generalized to the 
study of the dynamics of more interesting and non trivial examples but is not suitable for 
systems with collinear rods. 



1 These are structures for which the equations of statics do not suffice to determine all the internal forces. 
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FIG. 1. The four masses of magnitude m are connected with six rigid rods to form an equilateral 
triangle of side £y/3. The central mass is placed at the barycenter of the triangle. 


Our model consists of four equal masses in a triangular arrangement (the fourth mass 
is placed at the barycenter) connected by six rods as shown in Fig. [I] The system is 
constrained to move in the plane. Notice that the structure is rigid and remains so even 
after removing one of the rods. We take the following Lagrangian 

= IN| 2 Y ~ 0# “ > ( IL1 ) 

iev {i,j} 

where our notation makes use of the graph labels associated with our system according to 
the numbering shown in Fig. [lj Here the sets of vertices and edges are respectively given by 
V = {1,..., 4} and £ = {{1,2}, {1, 3}, {1,4}, {2, 3}, {2,4}, {3,4} j. We denote the length 
of the {i,j} edge as (see Fig. [l] for their values). The positions and velocities of the 
particles are represented by q* and v ?: respectively. Each of the configuration variables q{ij} 
is a Lagrange multiplier enforcing the condition that the length of the edge {i,j} is injy- 
These conditions are holonomic constraints (i.e. velocity independent) and, hence, can be 
introduced in the Lagrangian in this simple way without modifying the Newtonian equations 
of motion. The velocities associated with the qyijy are denoted as vyijy. 

We pause for a moment to compare L with other well known Lagrangians, in particular 
the one describing the free electromagnetic field given by 

L em ( A, Ho, V, Vo) = l - J ((V + VH 0 ) • (V + VH 0 ) - (V x A) ■ (V x A)) dx. (II.2) 

Here A and Ho are function^ on M 3 that denote the vector and scalar potentials respectively 
with the corresponding velocities being V and Vo- As we can see L and Lem share some 
features, for instance, the velocities associated with some of the dynamical variables {qyijy 
and H 0 ) do not appear. Although the introduction of the terms involving qujy in L may 
seem as an artificial complication they are, in fact, similar to the ones involving H 0 in Lem- 
We exploit this analogy in the present paper. 

The Euler-Lagrange equations derived from the Lagrangian L give 


= - Y ^{M'}(q*-Qj), iev, 

j~i 

q, - qj 2 = ^ J} , {i,j}e£, 


(II.3) 

( n -4) 


where the j ~ i notation means that the sum extends over all vertices j connected with the 
fixed vertex i. 

Notice that each term (q ? ; — q?) in the r.h.s of Eq. (II.3) can be interpreted as the 

force exerted on the particle i by the rod connecting it with particle j because the sum of 
these terms is the force acting on the ith particle. As there are no terms involving q{ij} the 
system is singular. It is easy to show the existence of families of solutions to these equations 


2 Their dependence on the spatial coordinates can be roughly understood as the presence of a continuous 
index in analogy with the index i £ V that labels the particles in our model. 
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that depend on a free function but describe the same physics. Consider for instance 


qi (t) = (o,o), 

(II.5a) 

q 2 (t) = (—f? sin cut, £ cos cot ), 

(II.5b) 

q 3 (f) = (—£sin(u;£ — 27t/3), £cos(u;t — 2 tt/3)) , 

(II.5c) 

qzi(t) = (—£sin(o;f + 27r/3), £cos(cut + 27r/3)) , 

(II.5d) 

5(1,2} (*) = 5(1,3} (*) = 5(1,4} (0 = fit) , 

(II. 6a) 

5(2,3} (t) = 5(2,4} (t) = 5(3,4} (t) = 7^(mu 2 /(t)) , 

(II.6b) 


where f(t) is an arbitrary function of time and oj a real parameter. It is straightforward 
to check that these functions satisfy Eqs. (II.3) and (II.4). Notice that the motion of each 
particle is perfectly determined, which implies that the force acting on each of them also is. 
However, the force exerted by each rod is not completely determined. This is exactly what 
happens in gauge theories: some magnitudes are arbitrary to a certain degree but there are 
“physical functions” of them (observables) that are not arbitrary. 

The constraints (II.4) imply that we cannot freely choose initial positions for all the 
particles in the system. They also imply that the velocities cannot be arbitrary either j^] 
A set of conditions that the velocities must satisfy can be obtained by differentiating Eq. 


(q. - Qj) • (qt - qj) = o, {*, j} e £. 


(II-7) 


At this point a simple procedure to obtain all the possible conditions on the configuration 
variables and velocities suggests itself: keep on differentiating and using, when possible, the 
equations of motion and the constraints already obtained to determine if new independent 
conditions appear. Although this method can actually be implemented in the present cas^] 
there are two reasons not to do so. The first is that it is somehow difficult to End out when 
the procedure stops. The second is that we are interested in the Hamiltonian formulation 
(necessary, for example, to attempt the canonical quantization of our system a la Dirac). 


III. THE FOUR PARTICLE MODEL: THE HAMILTONIAN PICTURE 


We obtain now the Hamiltonian formulation for our system by using a method inspired in 
the GNH algorithm [2j. The first step to get the Hamiltonian formulation for a mechanical 
model from its Lagrangian is to define the canonical momenta and write the generalized 
velocities in terms of them. Singular systems are identified, in practice, as those for which 
this is not possible. One might hastily conclude that the Hamiltonian formulation does not 
exist in this case; however, a quick look at the equations of motion suggests a possible way 
out: the fact that the positions of the particles and their velocities are subject to constraints 
such as (II.4) or (II.7) could mean that the proper space to define the Hamiltonian dynamics 
is not the “full phase space” T consisting of all the generalized positions q*, q{ t . 3 } and 


3 The velocity field for a rigid solid has a very specific form that determines the actual freedom to choose 
the initial velocities of the particles. 

4 The Lagrangian symplectic approach [3] would provide the rigorous geometric implementation of this 
procedure. 
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momenta p,, puj\ but rather an appropriate subset of it. The definition of the canonical 
momenta p = dL/dv from the Lagrangian L given by Eq. (II. 1), can be interpreted as a 
map (known in the technical literature [H] as the fiber derivative) FL : M 28 —> T defined by 


(q», Q{iJ}, v {ij}) ^ (qo dL/dvi, dL/dv {iJ} ) = (q i? mv*, 0). (111.1) 

As the momenta P{ijy are zero, the image under FL of a curve (q&(£), Qk{t), Q{i,j}(t)) 

must be contained in the so called primary constraint submanifold of T given by 


M 0 := {(q», q{i,j},Pi,P{ij}) 6 K 28 : C ( ^ j} := p {iJ} = 0 , {i,j} G 8} . (III.2) 

It is, hence, natural to look for a Hamiltonian description defined only on Mo or an appro¬ 
priate subset of it. Mo can be viewed as M 22 with coordinates (q,;, ?{m}, P?;)- 

The evolution of the system in Hamiltonian form is given by the integral curves parame¬ 
terized by time (q. t (t), q{ij}(t), p i(t),p {iij} (t)) of a vector field X = (X q ., X q{ . . } , X p ., X P{i< , } ). 
These are given by the first order differential equations 


q* Xqi ) Q{i,j} ^9{i,3'J 5 P i -^-q iip{i,j} ' (HI.3) 

The vector field X is obtained by using a construction that involves the Hamiltonian F[ and 
an important geometric object: the symplectic form fl. 

The usual prescription to obtain a Hamiltonian (based on performing a Legendre trans¬ 
form) cannot be applied to this example but a simple extension of it can be used on Mo- 
In the present case this amounts to ignoring the P{i,j}f{i,j} terms to get 

H{qi,q {itj} ,Pi) = ll pi H 2 + \ 9{ij}(||qi ~ Qil| 2 “ ■ ( IIL4 ) 

iev {i,j} e£ 

Although this function is dehned in principle only on Mo it can be trivially extended to T 
as H(qi,q {iij} ,pi,p {id y) = H(q h q {iJ} , p*). 

The inhomogeneous linear equation that would determine X if the system were non¬ 
singular 

nX-VH = 0 , (III. 5 ) 

where the 28 entries of f2X (dehned by the linear action of 81 on X) are real functions 
on the phase space and V FI denotes the 28-dimensional gradient of F[. The geometric 
structure of T is such that a specific symplectic form 81 with the required properties (non¬ 
degeneracy among them) can always be built —hence the word canonical. A famous theorem 
by Darboux [3] proves that it is always possible to find a coordinate system covering a large 
enough part of T (an open set) where 81 can be written in matrix form as 


( 0i4xl4 
V 14x14 


— Zl4xl4 
0l4xl4 


(III.6) 


This is actually the reason why elementary treatments avoid discussing the determination 
of the vector held X from the Hamiltonian through the solution of (III.5): when f2 takes 
the previous form the equations for the integral curves of X are the textbook Hamilton 


5 This form is enough for our needs in this paper. A more rigorous way to write this equation requires the 
use of differential forms and related concepts. 
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equations. The Hamiltonian treatment of singular systems requires in an unavoidable way 
the explicit consideration of the symplectic structure. This is why we mention it here. 

Given that the dynamics in our example must be confined to Ado (P{i,j} — 0), we can 
work as if this was the full phase space, in particular, try to find Hamiltonian vector fields 
X= (X,0) = (X q . , X q{ . . , X p ., 0), i.e. tangent to Ad 0 and defined only there. Acting on 
these vectors Eq. {[n becomes 


(uX-VH)\ Mo = 0, 
where u> is the 22 x 22 degenerate matrix 


(III-7) 


( 08x8 

08x6 

~hx8 

W = 1 06x8 

06x6 

06x8 

\hx8 

08x6 

08x8 


(III.S) 


The action of u> on X that we need in order to solve Eq. (III.7) is 

= (-X Pi ,0 6 ,X qi ). 

The gradient on Ad 0 is V := (<9 qi , <9 Pi ), hence, 


vi7= q* - Hi)’ ^(llq* - ^II 2 - • 


(III.9) 


(III.10) 




The equations (III.7) constitute a linear, inhomogeneous system so, generically, some condi¬ 
tion must be satisfied by the inhomogeneous term for the system to be solvable. It may also 
happen that only part of the unknowns (the components of X in this example) are fixed 
after solving it. In the present case we easily obtain 


X qi = - 

m 


X p = 

P i 


- qj), 




o = ||q* - qj| 


'{*j} 


= : C 


(i) 

{hi} 


i e V 
i e V 

{i,j} e £. 


(HI-11) 

(111.12) 

(111.13) 


As we can see Eq. (III.7) cannot be solved on the whole of Ado but only in the part of it 
where condition (III.13) is satisfied. Let us call this subset Adi, i.e. 


Ml := {(q®, e Ado : C^ j} = 0 ,{i,j}e£}. 


(III.14) 


The components of X must have the form given by Eqs. (III.11) and (III.12). Notice 
that the components remain arbitrary at this stage. Now, if the vector held X , 

with the form just obtained, were tangent to Adi we would have succeeded in Ending an 
appropriate submanifold of the phase space where we can define the Hamiltonian dynamics 
of our singular system. This can be shown by checking if the directional derivative 


of along X vanishes. In this case we have 


= Aq< - q. 


'{ijy 


m 


(p* - Pi) = : c {i,i }» 


(III.15) 
















( 2 ) 

which implies that X is only tangent to M\ at the points satisfying Cj/j = 0. These define 
the new submanifold 


M 2 ■= {(qn q{ijy,Pi) e Mi : = 0 , {z, j} G £} . 


(III.16) 


(o) ' . . 

The conditions Tj = 0 are necessarily independent of (III. 13) as they involve the momenta 
p, . It is interesting at this point to pause for a moment to understand their meaning. First 
of all, as the system is contained in the plane, it is obvious that 


(q* - qj) • (p i ~ P j) = o yy Pi - P j = u {itj} R(c[i - qj) , {i,j} g £, 


(III.17) 


where R is a counterclockwise rotation of n/2 and oo{ij} are real coefficients. If we select three 
of the particles in the system (say, 1, 2 and 3) and add the expressions given in Eq. (III. 17) 
for {i,j} = {1,2}, {2,3} and {1,3} we immediately see that U{ q 2 } = <^{2,3} = <^{ 1 , 3 } = : mu 
with u G M. Considering the remaining triangles in the graph associated with the system 
we get (npjj = mu for every {i,j} G £, that is, the constraints Cjy(, = 0 are equivalent to 
the existence of a real parameter u such that 


P* - Pj = muR(qi - qj) , {i,j} G £ . (III.18) 

As in the present case the velocities are just the momenta divided by m these last conditions 
are equivalent to saying that the velocities correspond to those of the particles of a rigid 
body (with an angular velocity given by u). 

The way to proceed is obvious now, compute 

VxC{^ } =^||Pi-p j || 2 -^g{i,fc}(q~q i Hqi-q fc )+^g{yfe}(qi-qiHqi-qfc)= :C {S} , (III.19) 

k^i k~j 


and check if the conditions C® — 0 for {i,j} G £ provide additional constraints. On the 
submanifold Ad 2 these conditions can be written in the form 


2 mu 


<2 

2 


e- 2 


- E 

{k,l}&£ 


Af, 


{k,l} 


{hj} 


Q{k,i} = 0 , {i,j}e£ 


(III.20) 


where, by using the geometry of the system (see Fig. 1), the matrix M = (M^ J .j fc ’^) can be 
seen to be 

/ 4 -1 -1 3 
-1 4-13 

-1-1 4 0 

3 3 0 12 


M = 


3 

0 

3 

3 


°\ 

3 


3 
3 

3 0 3 3 12 3 

VO 3 3 3 3 12/ 


(III.21) 


The entries of M are labelled in the order {1,2}, {1,3}, {1,4}, {2,3}, {2,4} and {3,4}. 
The rank of M is 5 and its kernel is spanned by the vector u T = (—3, —3, —3,1,1,1). It 
is straightforward to see that the conditions given by ( III.20[ ) are a compatible system of 
equations for the qujy where we can solve for any Eve of them in terms of the remaining 
one. The solutions of (III.20) can be parameterized, for example, in the form 


A mu 2 

<7{1,2} = <?{1,3} — 9(1,4} = A, 9(2,3} = 9(2,4} = 9(3,4} = H- 


A G 


(III.22) 


3 
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At this point we have the submanifold 


M 3 := {(q®, q{i,j},Pi) e M 2 : Cf\ } = 0, {i, j} e £} 


(III.23) 

(3) 


We need to check again if the vector field X is tangent to Xi 3 by requiring VxCl /4 = 0, 


i.e. 


o-E x nw ( qi - q i)' ( q?; - q*) - m E x nm (q*. - %•) • (q* - qfe) 


k~i 


k~j 


3 ~ Pi)-(q* - qfc) - 3 - Pj)-(qj - q*) (m.24) 


k~i 


k~j 


E - qj)-(Pi - P k) - E 9{i.fc}( q i - q?')'( p i - Pa 


k^i 


k~j 


By using again the geometry of the system and (jlll. 18|) we can write (|III.24|) in the form 

(III.25) 


E M {ij } + 2 ^ u E %/ 0 > e E 

{k,l}e£ {k,l}e£ 

where the matrix N = (TW^’^) is 

/ 0 -1 1 1-1 0 \ 

1 0-1-1 0 1 

-11001-1 
-1100-33 
1 0-1 3 0-3 

V 0 -1 1 -3 3 0 / 


N = 


(III.26) 


This is a linear inhomogeneous system of equations for the X q{ . j} components of the Hamil¬ 
tonian vector field X. No new constraints appear now as compatibility conditions because 
u T N = 0. We can then solve for the functions X q{ . to finally get 


1,2} — -^9{l,3} _ ^<?{1,4} 


Y _ Y _ Y 

^-<?{2,3} ^ V 9{2,4} ^^{3,4} 


(III.27) 


where we have used Eq. (III.22) and 5 is an arbitrary real function on Ai 3 . The algorithm 
stops here because there are no more conditions on the canonical variables and we have been 
able to solve for the components of the vector field X satisfying the tangency conditions 
and the basic equation (III.7). The final submanifold given by the algorithm is 

>(2) 


{(«,«(«>, Pi,«M»)er : C(S ) = C (.'i}= C S>= C S} = 0 '( i ’j} 6 £} (HI-28) 

and the A 43 -tangent vector field is given by (III.11), (III.12), (III.27). 

We check now that we get the dynamics described by the original equations of motion 


7(1) 


(3) 


(II.3) and (II.4). Indeed the equations for the integral curves of the Hamiltonian vector held 
X are 


Hi = 
P* = 

9(1,2} 

9(2,3} 


/m, 

i G V, 

(III.29) 


i e V, 

(III.30) 

9(1,3} = 9(1,4} = “ , 


(III.31) 

9(3,4} = 9(2,4} = ——/3 , 


(III.32) 
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with initial conditions satisfying the constraints 


C 


(i) _ r m 


~ Qu\ ~ ~ ^ ^ £ ' 


(III.33) 


= 0 is 


As we can see equations (III.29) and (III.30) imply ( |II.3 ). The constraint Cj /; j 
equivalent to (II.4) and the remaining constraints are necessary for the consistency of the 
dynamics (in particular to choose good initial data). As we can take S to be an arbitrary 
function on XI 3 and the initial data for the q{ij\ must satisfy (III.22), the solutions to 
(III.31) and (III.32) must have precisely the form given by Eqs. (II. 6 a), (II. 6 b). 

The presence of the arbitrary function E in the Hamiltonian vector field X is directly 
related to the gauge symmetry of our system. Suppose that we pick a point P 0 on the 
submanifold M. 3 , make several different choices of S, compute the integral curves of the 
resulting X starting from P 0 a t fo an d take the points of these curves corresponding to 
the same later value of the time parameter t > t$. From a physical point of view these 
configurations should be considered as equivalent (they certainly are, both regarding the 
positions of the particles and the forces acting on them at each instant of time). This leads 
us to the definition of gauge orbits as constituted by all the points in phase space reachable 
from allowed initial data after a certain fixed time by making any possible choice of the 
arbitrary part of the Hamiltonian vector field defining the dynamic^} In order to avoid the 
redundant description of equivalent physical configurations two options are available: gauge 
fixing and the introduction of the reduced phase space. We briefly describe them in turn. 

A popular way to select the arbitrary components of the Hamiltonian vector field X 
(encoded in S) is through gauge fixing. In our example this amounts to selecting the value 
of the force exerted by one of the rods. This can be realized physically by substituting one 
of them by a spring of fixed rest length (or even removing one rod). By introducing, for 
instance, the additional gauge fixing condition Q := ^{ 1 , 2 } = 0 we build a submanifold M. G 
of Ad 3 and fix S by demanding V xG — 0. This immediately gives S = 0 and removes the 
arbitrariness in the evolution. 

The reduced phase space is the abstract space of gauge orbits endowed with the appro¬ 
priate geometric structures (in particular a symplectic form and an appropriate restriction 
of the Hamiltonian vector fields defining the dynamics [3]). In the present case -but not in 
generic gauge theories such as electromagnetism, Yang-Mills or general relativity- it can be 
obtained by relying on the original idea by Lagrange to avoid constraint forces by working 
with appropriate “generalized coordinates” and writing the Lagrangian in terms of the ki¬ 
netic and potential energy. Here the appropriate coordinates are the position of the center 
of mass ( x,y ) € M 2 and an angle 6 , i.e. a point on the unit circle S 1 . The Lagrangian is 

3 


Lr(x, y, 9, v x , v y , v g ) = 2 m{v x + v y ) + -mv 9 . 
The system in this form is not singular, the Hamiltonian is 

H R (x,y,9,p x ,p y ,p g ) = ^-(pl + p 2 y ) + -^p 2 e 


giving the unique Hamiltonian vector field 


X R :=(X x ,X y ,Xg 


Xp x 5 Xp y , X pg 


= fPx_ Ihj_ Pl_ 
\4m’ 4m’ 3 m 


, 0 , 0,0 


(III.34) 


(III.35) 


(III.36) 


6 Gauge orbits can be characterized also geometrically by considering the degenerate directions of the 
“pulled back” symplectic form [3]. 
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It is important to mention that the reduced phase space, whose points are of the form 
(x,y,9,p x ,p y ,pg), is non-trivial as a manifold. Indeed it has the form T# = (M 2 x S 1 ) x M 3 
and, hence, is not isomorphic to a Euclidean space. 


IV. COMMENTS 


As we have shown it is possible to define simple mechanical models that behave as gauge 
systems in a non trivial way (some trivial examples can be found, for instance, in [5]). The 
main ideas are to implement the indeterminacy characteristic of hyperstatic structures in 
a dynamical setting and introduce constraint forces as explicit dynamical variables. The 
equations of motion for simple models consisting of point particles connected by ideal rods 
mimic the most important features of gauge theories. In this sense they provide a useful 
finite dimensional analogue of gauge held theories and help as good pedagogical models 
to discuss other important issues such as quantization (both in the Dirac approach and 
by using path integral methods). It is important to mention, nonetheless, that a concrete 
implementation of these mechanical systems would not be subject to any indeterminacy in 
the individual forces exerted by the rods as a consequence of their elastic properties. In this 
sense the gauge behavior that we have discussed is a feature of the equations of motion in 
the simplified setting where the elastic properties of the rods are neglected. As far as we 
can see there is no analogue of this phenomenon in the standard gauge held theories. 

We have discussed the obtention of the Hamiltonian description of the model by using a 
geometric approach. An interesting exercise is to derive the same results by following the 
standard method proposed by Dirac and based on the use of Poisson brackets. In our opinion 
the GNH approach that we have followed is both conceptually clean and easier to use. Had 
we not stopped to discuss the meaning of the conditions that we have been obtaining, the 
computation of the Hamiltonian vector held and the constraints for our model could have 
been written in one page. The description of the submanifold in the full phase space where 
the dynamics takes place is very economical: it is seen as an algebraic manifold defined by 
the vanishing of simple polynomials in the canonical variables. This means, in particular, 
that it is defined globally in a coordinate independent way. The key concept in the obtention 
of the relevant manifold where the dynamics is defined and the Hamiltonian vector held is 
tangency. This is both conceptually simple and easy to implement in practice. 

An interesting problem for the reader would be to consider the same system after removing 
one rod. In this case no gauge invariance remains despite the fact that the Lagrangian is 
still singular. This example can help in further clarifying the relationship between singular 
Lagrangians and the presence of gauge symmetries. 

Finally, though it may sound trivial in a sense, we would like to point out that we have 
adapted our notation to the graph naturally associated with our system. This approach 
may be useful to study the dynamics of more complicated models because the basic form 
of the constraints that we have obtained should generalize readily. Some of them may even 
be interesting as they would provide a novel way to study the rigidity of frames (a held 
where there are still open problems) and the dynamics of complex structures such as flexible 
polyhedra. 
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